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The transverse momentum (pr) spectra and ratios of identified charged hadrons (vr*, K^, p, 
p) produced in ^/sj^ = 200 GeV Au+Au and d+Au collisions are reported in five different cen- 
trality classes for each collision species. The measurements of pions and protons are reported up 
to pt = 6 GeV/c (5 GeV/c), and the measurements of kaons are reported up to pt = 4 GeV/c 
(3.5 GeV/c) in Au+Au (d+Au) collisions. In the intermediate pr region, between 2-5 GeV/c, a 
significant enhancement of baryon to meson ratios compared to those measured in p+p collisions 
is observed. This enhancement is present in both Au+Au and d+Au collisions, and increases as 
the collisions become more central. We compare a class of peripheral Au+Au collisions with a 
class central d+Au collisions which have a comparable number of participating nucleons and binary 
nucleon-nucleon collisions. The pr dependent particle ratios for these classes display a remarkable 
similarity, which is then discussed. 

PACS numbers: 25.75.Dw, 25.75.Ld 



I. INTRODUCTION 

Measurements of identified particles in Au+Au colli- 
sions allow the study of particle production mechanisms 
in a hot and dense nuclear medium and probe the prop- 
erties of the quark-gluon plasma (QGP) [H-Q. In d+Au 
collisions, these measurements allow the study of cold nu- 
clear matter effects onparticle production, such as the 
Cronin enhancement [J Q , nuclear shadowing Q , and 
gluon saturation [1]. These cold nuclear matter effects 
are present in Au+Au collisions as well, and the study 
of d+Au collisions allows us to determine these effects 
directly and to disentangle them from the effects of the 
hot and dense nuclear medium. 

One of the most intriguing discoveries in the early days 
of the research program at the Relativistic Heavy Ion Col- 
lider was the significantly enhanced baryon production 
relative to meson production at intermediate transverse 
momentum 2 GeV/c < pr < 5 GeV/c, as evidenced in 
the large baryon to meson ratios and the significant dif- 
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ferences in the particle suppression patterns measured by 
the nuclear modification factors [9l-[llj . Several classes of 
models were introduced to explain these differences based 
on different physical phenomena, such as strong radial 
flow [l^j E^l' baryon junctions fl3 . [TBj . and hadroniza- 
tion through recombination [16l42l| . 

Additionally, the recombination model has been em- 
ployed to explain the baryon vs. meson difference in the 
Cronin enhancement observed in p+ A collisions [12, ■ 
Traditional explanations of the Cronin enhancement in- 
volve the multiple soft scatterings in the initial state 
prior to the hard scattering and subsequent fragmen- 
tation of the hard scattered parton [24i|. This process 
can naturally explain the deficit of particle production 
at low pt and enhancement at intermediate pt, but 
does not account for the particle species dependence at 
Relativistic- Heavy-Ion-CoUider energies [l^. Conversely, 
in the recombination model, the observed enhancement 
is attributed to final state effects, i.e. the recombina- 
tion of soft partons from the nuclear medium with hard 
scattered partons in a jet. For this reason d+Au col- 
lisions represent an excellent testing ground for the re- 
combination model, since hot nuclear matter effects, such 
as the collective expansion of the medium, are not ex- 
pected. However, recent results in d+Au collisions at 
200 GeV di] and p+Pb collisions at 5.02 TeV [13^ 
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suggest some collective expansion effects may be present 
in the most central events. 

Measurements of strange particles, such as charged 
kaons, have also been an interesting subject in heavy ion 
collisions. An enhancement of strangeness production 
relative to that in p+p collisions has been observed at 
various colliding energies (soj . This strangeness enhance- 
ment is a possible signature of deconfinement, thermal- 
ization, and flavor equilibration (3ll.[33|. In this scenario, 
strangeness production is dominated by thermal gluon 
fusion. The measurement of charged kaons in a broad 
Pt range and in different centrality classes is a signifi- 
cant tool to further understand the thermalization of the 
system and the mechanism of strangeness production. 

To address the particle production in both hot and 
cold QCD matter, a systematic study of identified parti- 
cles over a broad pt range with a wide selection of cen- 
tralities in both Au-f An and d-t-Au collisions is required. 
In this paper, the spectra, particle ratios, and nuclear 
modification factors previously reported by PHENIX in 
Au-f Au [l0| and d+Au collisions [1^ are revisited, ex- 
tending the Pt reach of previous measurements and sig- 
nificantly improving the statistical precision. In Sec- 
tion |lT] we discuss the experimental apparatus and the 
detector subsystems used in this analysis; in Section Hill 
we discuss the analysis method, including event and track 
selection, and particle identification; in Section IIVI we 
discuss the results; and in Section El we summarize our 
findings. 
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FIG. 1: (Color online) The PHENIX detector configuration 
for the 2007 data taking period. The 2008 configuration is 
very similar. 



II. EXPERIMENTAL SETUP 

The PHENIX experiment is a large, general purpose 
detector with a wide variety of detector subsystems ide- 
ally suited to the study of nuclear matter in conditions 
of extreme temperature and density. PHENIX is com- 
posed of global event property detectors, forward and 
backward rapidity arms (North and South) dedicated to 
muon measurements, and two central arm spectrometers 
(East and West) at midrapidity covering pseudorapidity 
region of |?7| < 0.35 for measurements of photons, elec- 
trons, and charged hadrons. Detailed descriptions of the 
various detector subsystems can be found in [s^ . 

Figure m gives a schematic diagram of the PHENIX de- 
tector for the 2007 configuration; the 2008 configuration 
is very similar. The upper panel shows the central spec- 
trometer arms, viewed along the beam line facing North. 
The lower panel shows the two forward rapidity muon 
arms (North and South) and the global detectors. 

The analysis presented herein makes use of the beam- 
beam counters, the tracking system (drift chamber 
and two layers of pad chambers), the electromagnetic 
calorimeter, and the time-of-flight detector in the West 
arm. 



A. Detector subsystems 

The beam-beam counters (BBC) are used for the 
minimum-bias trigger, the centrality definition, the de- 
termination of collision vertex along the beam axis (the 
2:-vertex), and the event start time. The BBCs [b^. IsBj 
are located at ± 144 cm from the nominal interaction 
point. They cover 27r in azimuth and a pseudorapid- 
ity range of 3.0 < |?7| < 3.9. Each BBC is an array of 
64 identical hexagonal detector elements, with the beam 
pipe passing through the center of the array. Each el- 
ement is a quartz Cerenkov radiation counter, and the 
radiator and photomultiplier tube are constructed as a 
single piece. 

The PHENIX tracking system is optimized for the high 
multiplicity environment of ultra-relativistic heavy ion 
collisions. It comprises drift chambers (DC) [ssj and 
pad chambers (PCI, PC2, and PC3) [37l]. This analy- 
sis makes use of the first layer of pad chambers (PCI) 
and the third layer (PC3). The DCs have an active vol- 
ume in the radial range 2.02 m < r < 2.46 m. The PCI 
is mated directly to the back of the DC frame in each 
arm at a radial distance of 2.49 m. The PC3 is located 
in each arm at a radial distance of 4.98 m. 

The primary particle identification detector used in 
this analysis is the time-of-flight detector in the West 



5 



half of the central arm spectrometer (TOFW). The 
TOFW [M Hi is located at a radial distance of 4.81 m 
from the interaction point, having pseudorapidity cover- 
age of I77I < 0.35 and azimuthal coverage of 22° in two 
separate sections. The individual elements are multi- 
gap resistive plate chambers (MRPCs). Each MRPC 
has six 230 ^m gas gaps separated by five 550 /im thick 
glass plates. On each side of the outermost glass plates 
(1.1 mm thick) are carbon tape electrodes held at +7 kV 
on one side and -7 kV on the other side for a total bias 
voltage of 14 kV. When charged particles traverse the 
detector, the gas between the plates is ionized and the 
image charge is collected at each side of the chamber on 
four copper readout strips. Each strip has dimensions 37 
X 2.8 cm^ with a separation of 0.3 cm between them. The 
strips are oriented lengthwise along the azimuthal direc- 
tion. Each strip is read out from both top and bottom 
so that the time difference between them can be used to 
determine the hit position along the length of the strip 
with resolution of order 1 cm. The TOFW system is 
composed of a total of 128 MRPCs, 512 strips, and 1024 
readouts. The total timing resolution, which includes the 
uncertainty in the start time from the BBC, is 84 ps in 
Au-I-Au collisions In d-l-Au collisions it is 95 ps, 

where the slightly poorer resolution is due to the lower 
resolution of the start time determination from the BBC. 
This is due to the lower multiplicity in d+Au collisions 
and, for the same reason, the z-vertex resolution is also 
poorer in d+Au collisions. 



III. ANALYSIS METHOD 
A. Event and Track Selection 

This paper presents an analysis of Au-|-Au collisions at 
y/s^ = 200 GeV, collected in 2007, and d-fAu collisions 
at y/sNN = 200 GeV, collected in 2008. For each set we 
select events that pass the minimum-bias trigger, which 
is defined as a coincidence between the North and South 
BBCs. In Au+Au collisions, this trigger requires two or 
more photomultiplier tubes firing in each BBC and mea- 
sures 92±3% of the total inelastic cross section; in d+Au 
collisions, it requires one or more photomultiplier tubes 
firing in each BBC and measures 88±4% of the total in- 
elastic cross section. We have an additional requirement 
that the collision vertex is within \z\ < 30 cm of the 
nominal origin of the coordinate system. 

Centrality selection is performed with the BBCs using 
the Glauber Monte Carlo procedure described in 40] , in 
which the charge in each BBC detector is assumed to 
be proportional to the number of participating nucleons 
A^part traveling towards it. For the Au+Au system the 
north and south BBC distributions are summed, but for 
the d-l-Au system, only the south (Au-going) side is used. 
The BBC charge is assumed to follow a negative binomial 
distribution (NBD) with a mean of Afpart and the remain- 
ing NBD parameters determined from a minimization 



of the combined Glauber-(~NBD calculation with respect 
to the data. The BBC distributions are divided into equal 
probability bins, and the corresponding Glauber distri- 
butions are used to calculate A^part as well as the number 
of binary nucleon-nucleon collisions A"coii, as shown in 
Table H 



TABLE I; Values of (Vcoii) and (Apart) for Au+Au and d+Au 
collisions from Glauber model simulations. 



Centrality 


(Vcoll) 


(Apart) 


Au+Au 






0-10% 


960.2 + 96.1 


325.8 + 3.8 


10-20% 


609.5 + 59.8 


236.1 + 5.5 


20-40% 


300.8 ± 29.6 


141.5 + 5.8 


40-60% 


94.2 ± 12.0 


61.6 + 5.1 


60-92% 


14.8 ± 3.0 


14.7 + 2.9 


d+Au 






0-20% 


15.1 ± 1.0 


15.3 + 0.8 


20-40% 


10.2 ± 0.7 


11.1 + 0.6 


40-60% 


6.6 + 0.4 


7.8 + 0.4 


60-88% 


3.1 ± 0.2 


4.3 ± 0.2 


0-100% 


7.6 ± 0.4 


8.5 ± 0.4 



Charged track reconstruction in the DC is based on 
a combinatorial Hough transform, which gives the angle 
in the main bend plane {r-(j)) and thus px- The PCI 
is used to determine the hit position in the longitudi- 
nal (z) direction. Only tracks with valid information in 
both the DC and PCI are used in this analysis. Tracks 
in DC/PCl are projected to the outer detectors, such as 
PC3 and TOFW, and matched to hits in those detectors 
with the minimum distance between the projection and 
the hit position. The distribution of differences between 
hits and projections is approximately Gaussian, with an 
additional underlying background caused by random as- 
sociations. Only tracks with a difference of less than 
two standard deviations in both the azimuthal and lon- 
gitudinal directions in both the PC3 and the TOFW are 
selected, so as to minimize background contamination. 
In the Au+Au data for pr > 5.0 GeV/c, an additional 
background isolation cut is applied. For these tracks we 
require E/px > 0.2, where E is the energy deposited in 
the electromagnetic calorimeter. This cut removes low pT 
particles that are falsely reconstructed as high pT tracks. 



B. Particle Identification 

Charged particle identification (PID) is performed by 
simultaneous measurement of the momentum, time-of- 
flight, and path- length. These quantities are used to de- 
termine the mass of the candidate track based on the 
following relationship: 
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where m is the mass, p is the momentum, c is the 
speed of hght, t is the time of flight, and L is the path- 
length. The distributions are approximately Gaus- 
sian. The standard deviation of the distribution can 
be parametrized as a function of momentum as follows: 

-l^ - f (4-V) + |f(4m^(l + ^)) 

+ (2) 

where tTm2 is the standard deviation of the rri^ distribu- 
tion, m denotes the physical mass of the particle and thus 
the square root of the centroid of the rn^ distribution, CTq, 
is the angular resolution of the drift chamber, is the 
multiple scattering term, at is the total timing resolu- 
tion, and Ki is the magnetic field integral constant. The 
magnetic field integral constant depends on the magnetic 
field configuration. The PHENIX magnet system for the 
central arms comprises two coils. The two coils can be 
run together, opposed, or with the inner coil off. During 
the 2007 Au-t-Au data taking, the coils were run opposed, 
while for the 2008 d-t-Au data taking, the coils were run 
together. Running the coils opposed produces zero mag- 
netic field in the region between the beam-pipe and the 
inner coil of the magnet. This is needed for the analy- 
sis of the dielectron continuum using the hadron blind 
detector, which is the innermost detector during this op- 
erational period and can be seen (labeled as HBD) in 
Fig. [1] The PID parameters are presented in Table HIl 



TABLE II: Parameters for the PID function defined in Equa- 
tion 



Parameter 




2007 Au+Au 


2008 d+Au 




(mrad) 


0.896 


1.050 




(mrad-GeV/c) 


0.992 


1.000 




(ps) 


0.084 


0.095 


K, 


(mrad-GeV/c) 


75.0 


104.0 



To select candidate tracks of a particular particle 
species, the is required to be within two standard de- 
viations of the mean for the selected particles species and 
outside two standard deviations of the mean for the other 
particle species. Below the regions where the PID bands 
intersect, the PID contamination is negligible. Above it, 
contamination is from one side of the tail of the distri- 
bution beyond 2 standard deviations. For a Gaussian 
distribution this is 2.25%. Therefore the contamination 
can be estimated based on the ratio of the respective 
yields multiplied by this value. The distributions 



may have slightly nonGaussian tails and therefore the 
PID contamination may be slightly higher. In any event, 
we estimate the PID purity to be better than 90% for 
all particle species at all pr, in all centrality classes, and 
in both collision systems. Figure [2] shows the vs the 
PT multiplied by the charge for the 2007 Au-f Au data; 
the 2008 d+Au data are very similar. The 2a PID bands 
are superimposed as solid black lines. The upper panel 
shows the entire distribution with all the track selec- 
tion cuts applied but none of the PID cuts. The lower 
panel shows the same distribution but with the PID cuts 
also applied. 




-6 -4 -2 2 4 
charge x (GeV/c) 



FIG. 2: (Color online) The particle identification method, m? 
vs charge x pr for 2007 Au+Au data; the 2008 d+Au data are 
very similar. The solid black lines indicate the two standard 
deviation PID bands used for the cuts. The upper panel shows 
the bands superimposed on the entire distribution, the 
lower panel shows the distribution after the cuts have been 
applied. 



C. Corrections to the Raw Data 

In order to obtain the true invariant yield, the raw 
spectrum needs to be corrected for a variety of factors. 
Various types of simulations are performed to determine 
these corrections. To correct for geometrical acceptance. 
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analysis cuts, particle interactions with detector mate- 
rials, and in-flight decays (for pions and kaons), we use 
single particle Monte Carlo (MC) simulations. For these 
simulations, single particles are generated using a random 
generator, with flat distributions in rapidity, azimuth, 
and pt- The random particles are then run through a 
GEANT simulation of PHENIX to determine the interac- 
tions of the single particles with the detector subsystems 
and support structures. Next, all the detector response 
information is fed through the usual PHENIX reconstruc- 
tion software to produce simulated tracks. Finally, these 
simulated tracks are analyzed in the exact same way as 
tracks from the real data in order to determine the cor- 
rections. The total correction factor, Fc{pt), is given by 
the following relation: 

p / N dNputput/dpT _ . ^ 

^C\PT) — / , — ^acccptancc^cfficicncv^cuts • 1<J/ 

dNin-put/dpT 

To correct for the detector occupancy effect, which is 
most signiflcant in the TOFW, we run embedding sim- 
ulations, where a track from single particle MC simula- 
tions is embedded into a real event, and the occupancy 
correction is determined from the relative efficiencies of 
reconstructing the single track in isolation and in the 
event. This correction is the largest in the most central 
Au-|-Au collisions where the multiplicity is the highest 
and therefore the occupancy effect is the strongest. In 
the most peripheral Au-|-Au collisions the multiplicity is 
low enough that there is essentially no effect. The same 
is true in d-l-Au collisions, where no correction is applied. 

For the d-|-Au system we apply a correction for the un- 
derlying event correlation that exists between produced 
particles measured in the central arms and particles at 
forward angles that satisfy the BBC interaction trig- 
ger [4l|. This correlation produces both a trigger bias, 
in which events satisfying the trigger are biased towards 
higher multiplicities, and a bin shift, in which nominally 
peripheral events are shifted to higher centrality bins, 
thereby depleting the more peripheral bins. We correct 
for these effects using the Glauber Monte Carlo combined 
with central particle yields measured in p+p collisions. 
Using this same framework, we also generate a correc- 
tion factor to convert the minimum-bias sample (0-88%) 
into one with zero-bias (0-100%). 

Table IIIII shows the centrality dependent corrections 
for each collision species. The occupancy is represented 
as an efficiency, while the bias factor is represented as a 
multiplicative correction. 

The proton and antiproton spectra are additionally 
corrected for the feed-down from weak decays of hyper- 
ons into protons. We use single particle MC simulations 
of the A baryon and apply the analysis cuts used for the 
protons to determine the percentage of A baryons that 
decay into protons that pass our proton selection cuts. 
This is used to determine the percentage of the total pro- 
ton sample that likely comes from hyperon decays, which 
is called the feed-down fraction. 



TABLE III: Summary of centrality dependent corrections. 



Au+Au Occupancy Dependent Efficiency 

0-10% 0.542 

10-20% 0.653 

20-40% 0.783 

40-60% 0.904 

60-92% 0.964 

d+Au Bias Factor Correction 

0-20% 0.94 

20-40% 1.00 

40-60% 1.03 

60-88% 1.03 

0-100% 0.89 



The feed-down fraction is only dependent on the K/p 
ratio and not explicitly on the A spectrum itself. For the 
Au-|-Au data we take the K/p ratio to be 0.89 and the 
K/p ratio to be 0.95, as measured in Au-|-Au collisions 

y/s^ =130 GeV 42]. These are very similar to (and 
well within the systematic uncertainties of) the values 
of 0.91 for A/p and 0.94 for A/p obtained for Au-f Au 
collisions at ^/s^ = 200 GeV by examining the dN/dy 
values for A and A and p and p ^Td\ . We note however 
that the A data reported in Ref. [4^ are inclusive while 
the data A reported in Ref. [i^l are corrected for feed- 
down from other hyperons. 

For the d-l-Au data, we take the A/p ratio to be 0.85 
and the A/p ratio to be 0.99. There are no published 
data on A/p and A/p ratios nor A and A spectra or yields 
in d-l-Au collisions that could be used to estimate these 
ratios. Instead, we use the {A+A)/{p+p) ratio measured 
in p+p collisions [3], the p/p ratio measured in d-l-Au 
collisions [1^, and an estimate of the A/ A ratio in d+Au 
collisions based on measurements in Au-f Au [45j and p+p 
collisions [46j . 

The spectral shape of the A (A) is assumed to follow 
the p (p) spectrum with niT scaling. We also take the 
ratio to be independent of centrality. In fact, a small 
centrality dependence can be seen when examining the 
integrated yields in Refs. |Tol] and [43*1, although the dif- 
ferent centralities are consistent within the systematic 
uncertainties. 

All measurements of A implicitly include the S'^, which 
decays electromagnetically with 100% branching ratio 
to the A and a photon. We do not correct for feed- 
down from the charged S states, nor for the 5 and 51 
multistrange baryon states. Because these corrections 
are smaller and have large uncertainties, they are in- 
cluded in the overall systematic uncertainty estimates. 
For the charged S states we are concerned only with the 
S+, which decays to a proton with a 0.5 branching ra- 
tio. The data from 200 GeV p+p collisions show a ra- 
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tio (S+ + E- + E+ +_S-)/(A +_A) of 0.50 ± 0.18(syst) 
Ref. [Hi. Assuming S/E « A/A and E-/E+ « 1, 
we can estimate E+/A to be 0.25. The ratios for 
(=0 E°)/i^A + A) and + S+)/(A + A) in p+p colli- 
sions both have a value of 0.065 with very large («100%) 
systematic uncertainties. Assuming S/S « A/ A one 
can estimate the (S*^ + S~)/A ratio to be roughly 0.13, 
which is consistent with the total feed-down correction of 
15% reported for A in Ref. [i^. The latter also includes 
feed-down from 17, though the contribution is small. 

The systematic uncertainty on the feed-down fraction 
is estimated to be 25%, which is due primarily to the un- 
certainty in the A/p ratio. We also note that, in the pr re- 
gion of interest to this paper, 2 GeV/c < pr < 5 GeV/c, 
the feed-down fraction is of order 10%, so a 25% change 
in the A/p ratio produces only a roughly 2.5% change in 
the proton spectrum. 

Figure [3] shows the feed-down fraction as a function 
of pt for protons as solid black lines and antiprotons as 
dotted red lines for the 2007 Au+Au data (left panel) 
and 2008 d-l-Au data (right panel). As mentioned above, 
the magnetic field configuration is different for 2007 and 
2008, which is the reason for the different feed-down frac- 
tions between the two data sets. 



2008 -P 
d+Au p 




(GeV/c) 



(GeV/c) 



acceptance, track selection, and PID, are predominantly 
type B. The normalization corrections for effects such as 
detector occupancy and efficiency are type C. The uncer- 
tainties are assessed for each cut by redoing the analysis 
with the cut varied and then determining the difference. 
The cuts are varied in exactly the same way for the anal- 
ysis of both the experimental and the simulated data. 
This process is repeated for all the analysis cuts and the 
differences are summed in quadrature to determine the 
final uncertainty. This is done for both the spectra and 
for the ratios. The uncertainties are examined in each 
centrality class for both Au-|-Au and d-|-Au collisions. In 
Au-f Au collisions, they are found to be quite similar in 
all centrality classes. In d-l-Au collisions they are found 
to be negligibly different for all centrality classes. When 
taking ratios of the various particles, all of the type C 
uncertainties and some of the type B uncertainties can- 
cel. For antiparticle to particle ratios, the uncertainty 
from acceptance does not cancel at all, and the uncer- 
tainty from track selection and PID mostly cancel. For 
other particle ratios, K/tt and p/tt, the uncertainty from 
acceptance and track selection mostly cancel, and the 
uncertainty from PID does not cancel at all. In this 
analysis, we find the remaining systematic uncertainty 
on each of the particle ratios is roughly 5% for all pr- 
For the nuclear modification factor Rcp, which com- 
pares two different centrality bins of the same particle 
species in the same collision system, all the type B sys- 
tematic uncertainties cancel almost completely. There is 
an uncertainty of about 2% based on small variations of 
the track matching and PID distributions as a function 
of centrality. The type C uncertainties are completely 
uncorrelated and added in quadrature. For the nuclear 
modification factors Raa and RdA, the p+p reference 
data were collected during a different operational period 
and using different detector subsystems, therefore none 
of the systematic uncertainties cancel. A summary of the 
type B systematic uncertainties for the spectra is given 
in Table irTl 



FIG. 3: (Color online) Feed-down fraction as a function of 
Pt for protons (sold black lines) and antiprotons (dotted red 
lines) for the 2007 Au-|-Au data (left panel) and the 2008 
d-|-Au data (right panel) . The shaded bands indicate the sys- 
tematic uncertainty of 25%. 



D. Estimation of Systematic Uncertainties 

The different types of systematic uncertainty are cat- 
egorized as follows: type A, point-to-point uncorrelated 
in Pt, where the points can move up or down indepen- 
dently of each other; type B, point-to-point correlated 
in Pt, where the points can move together changing the 
shape of the curve; and type C, an overall normalization 
uncertainty in which all points move up or down by the 
same factor. The corrections for analysis cuts, including 



TABLE IV: Summary of systematic uncertainties from ac- 
ceptance, track selection, and PID of invariant yield of each 
particle species. 





71 + 


TV 


K+ 


K- 


P 


P 


2007 Au-I-Au 














PT < 3 GeV/c 


9% 


9% 


11% 


11% 


10% 


10% 


PT > 3 GeV/c 


10% 


10% 


11% 


11% 


11% 


11% 


Pt > 5 GeV/c 


14% 


14% 






14% 


14% 


2008 d+Au 














PT < 3 GeV/c 


8% 


8% 


13% 


13% 


9% 


9% 


PT > 3 GeV/c 


9% 


9% 


13% 


13% 


11% 


11% 



The type C uncertainties are from the centrality depen- 
dent corrections for the spectra and from the uncertainty 
on the Glauber model calculations for the nuclear modifi- 
cation factors. The uncertainty on the occupancy correc- 
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tion for Au+Au collisions is roughly 10%, and therefore 
the uncertainty on the yield varies from 5% for the most 
central to less than 1% for the most peripheral. The 
uncertainties for the Glauber values for A^coii are much 
larger and therefore dominate the uncertainty in the nu- 
clear modification factors. The uncertainty for the bias 
factors for d+Au collisions varies from about 1% for the 
most central bin to about 5% for the most peripheral 
bin. The bias factors are determined in the same Glauber 
model analysis as the iVcou and iVpart values, and there- 
fore the uncertainties are correlated. The uncertainty on 
the ratio of bias factors and iVcoii values used to deter- 
mine the nuclear modification factors varies from about 
3% for the most central to about 8% for the most periph- 
eral. 



IV. RESULTS AND DISCUSSION 

A. Invariant Yield as a Function of Transverse 
Momentum pt 

The main result of this study is the measurement of the 
invariant yield of pions, kaons, and protons as a func- 
tion of Pt in different centrality classes. The central- 
ity classes studied in the Au-|-Au measurement are 0- 
10% (the most central), 10-20%, 20-40%, 40-60%, and 
60-92% (the most peripheral). For the d+Au measure- 
ment, the centrality classes are 0-20% (the most central), 
20-40%, 40-60%, 60-88% (the most peripheral), and 0- 
100%. From these quantities, all other observables are 
derived, such as particle ratios and nuclear modification 
factors. Figure m shows the invariant yields of positive pi- 
ons, positive kaons, and protons (upper left, middle, and 
right, respectively), and negative pions, negative kaons, 
and antiprotons (lower left, middle, and right, respec- 
tively). The yields are scaled by arbitrary factors indi- 
cated in the legend for the sake of clarity and to keep the 
collision species grouped together. 



B. Particle Ratios as a Function of Transverse 
Momentum 

One of the simpler classes of derived quantities is the 
antiparticle to particle ratio. In the present analysis 
those ratios are 7r~/7r+, / , and p/p, which are 
plotted as a function of pt in the upper, middle, and 
lower panels of Fig. [SI respectively. In each panel, the 
Au-I-Au data are on the left and the d-fAu data are 
on the right. Drawn as a visual aid are dashed black 
lines with value 1.0 for the pions, 0.93 for the kaons, and 
0.73 for the protons; the values for kaons and protons 
are chosen from the reported pt integrated values from 
Ref. [l3l . Shown as a reference are data from p+p colli- 
sions from Ref. [13] for pr < 3 GeV/c and from Ref. Q 
for Pt > i GeV/c. Remarkably, all the ratios are es- 
sentially independent of both pT and centrality. Based 



on simple arguments about isospin conservation and the 
basics of the parton distribution functions and fragmen- 
tation functions, one would expect each of the antiparti- 
cle to particle ratios to vary as a function of pt as dis- 
cussed in Ref. [1^ . Indeed, these ratios have a significant 
Pt dependence in p+p collisions at midrapidity, and the 
agreement with theory depends significantly on the frag- 
mentation functions used [1^. However, the pt range 
needed to observe the decrease in these ratios in p+p col- 
lisions is quite large. As seen in Fig.O the pt dependence 
of the ratios in p-\-p is small and may be consistent with 
the ratios in d-|-Au and Au-|-Au. 

Figures [6] and [7] show the kaon to pion ratios as a 
function oi pt {K^ /tt'^ on the top, K~ /tt^ on the bot- 
tom) in Au-I-Au and d-|-Au collisions, respectively. The 
ratios in Au-I-Au collisions show a significant increase 
with increasing pT and a small increase as the colli- 
sions become more central. The enhancement of the 
integrated K/n ratio in more central collisions is at- 
tributed to strangeness equilibration in various thermal 
models [13, [HI, which is reflected in the differential ratio. 
However, the differential ratio may include additional 
information about the differences in the fragmentation 
functions and/or the phase space distribution functions 
used in the recombination models. As discussed in a pre- 
vious PHENIX publication [52| , the strangeness enhance- 
ment present in the hot and dense nuclear medium has 
an effect on certain recombination models [s^ . These re- 
combination models involve the recombination of partons 
in dissimilar momentum space, meaning that a shower 
parton from a jet can recombine with a thermal parton in 
the medium. The thermal component of thermal+shower 
recombination is more dominant at higher pT for strange 
hadrons (like kaons) than it is for nonstrange hadrons 
(like pions), leading to an enhancement of the ratio that 
increases with pT- This increasing enhancement mani- 
fests as the ratio rising more quickly in Au-I-Au collisions 
compared with p-\-p collisions, which is seen in Figs. [B] 
and [71 At sufficiently high pT where the shower compo- 
nent begins to dominate for both strange and nonstrange 
particles, this ratio is expected to turn over and begin to 
decrease. However, this turnover point, if it exists, is 
beyond the pt reach available for kaons in this study. 

The K/ TT ratios in d-f Au collisions are essentially iden- 
tical for all centrality classes, which may indicate that the 
mechanism for strangeness production in d-|-Au collisions 
is the same for all centrality classes. However, we also 
note that the various d-l-Au centrality classes span a rel- 
atively small range of A^part ■ Therefore, if the strangeness 
enhancement is only weakly dependent on A^part , the vari- 
ation of A^part in the d+Ku centrality classes may not be 
large enough for an effect to be observed. 

Figures [H] and [S] show the proton to pion ratios as 
a function of pT {p/t^^ on the top, p/ir^ on the bot- 
tom) in Au-I-Au and c?-|-Au collisions, respectively. Note 
that for Fig. [3| the vertical scale is different. The ratios 
in central Au-I-Au collisions show a strong enhancement 
over the values in p+p collisions. This is conjectured to 
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FIG. 4: (Color online) Invariant yield of tt*, K^, and p and p as a function of pr in Au+Au and d+Au collisions. The yields 
are scaled by the arbitrary factors indicated in the legend, keeping collisions species grouped together. 
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FIG. 5: (Color online) Ratio of invariant yield of -k~ /n'^ (top), K~ /K^ (middle), and p/p (bottom) as a function of pt in 
Au+Au (left panels) and d+Au collisions (right panels) in each centrality bin. Dashed lines are drawn as a visual aid with 
values of 1.0 for 7r~/7r"'", 0.93 for K~ / K'^ , and 0.73 for p/p. These values are taken from Ref. [l^. Shown as a reference are 
data from p+p collisions from Ref. ;47i] for pt < 3 GeV/c and from Ref. [4§| for pt > 3 GeV/c. 
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FIG. 9: (Color online) Ratio of invariant yield of protons 
to positive pions (upper panel) and antiprotons to negative 
pions (lower panel) as a function of pr in d+Au collisions 
in the centrality bins marked in the legend. Data for p+p 
collisions [431 are shown as a reference. 
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be attributed to the parton recombination mechanism of 
hadronization, which gives rise to a significant enhance- 
ment of baryon yields relative to meson yields in heavy 
ion collisions [la, [l^ Hoj. The p/ir ratios in the other 
centralities in Au+Au collisions show a clear and consis- 
tent trend with decreasing enhancement as the collisions 
become more peripheral. In d-f Au collisions there is a 
similar trend. The p/tt ratio in the most central d-|-Au 
collisions appears consistent with the ratio in the most 
peripheral Au-|-Au collisions. Additionally, the p/ir ratio 
is enhanced over p+p collisions for each centrality class 
in d+Au collisions except for the most peripheral. 



C. Nuclear Modification Factors as a Function of 
Transverse Momentum 

To measure the modification of the spectrum of pro- 
duced particles in heavy ion collisions relative to the spec- 
trum in p+p collisions, nuclear modification factors are 
employed. The nuclear modification factor i?AA is de- 
fined as the yield in Au-|-Au collisions divided by the 
yield in p-\-p collisions, normalized by the number of bi- 
nary nucleon-|-nucleon collisions iVcou as determined from 
the Glauber model. The nuclear modification factor Rcp 
is defined as the yield in central Au-f-Au collisions divided 
by the yield in peripheral Au-f-Au collisions, normalized 
to the respective numbers of binary nucleon-|-nucleon col- 
lisions. These can be expressed mathematically as: 



-Raa = 



{dN/dprY 



N^^^^idN/dpTy+P 



Rcp — 



{dN/dpTY^'^''-''' jyp^'^j^p''^'-" 

(dN / dpTY'^"^^'^^"'^ N'^'^"*^"-^ 



(4) 



(5) 



Figure M shows Rcp for 0~10%/40-60% (left panel) 
and 0-10%/60-92% (right panel) as a function of pt 
for charge averaged pions, kaons, and protons. Both pi- 
ons and kaons exhibit a suppression pattern at all val- 
ues of Pt- The kaons exhibit less suppression than the 
pions, indicating the additional role of strangeness en- 
hancement in the particle production mechanism. The 
observed enhancement of kaons relative to pions appears 
to be lower for the 0-10%/40~60% as compared to the 0- 
10%/60-92%, suggesting a centrality dependence of the 
strangeness enhancement, as seen in the K/n ratios dis- 
cussed above. The protons on the other hand exhibit 
quite different behavior, rising to a value very close to 
unity, indicating nonsuppression, around 2-3 GeV/c in 
PT- At higher values of pT the proton Rcp falls off slowly, 
beginning to approach the pion Rc p at the highest values 
of Pt available. The proton Rcp shown here is consis- 
tent within the systematic uncertainties with the proton 
Rcp reported by STAR 

Figure [TT] shows Raa as a function of pt in different 
centrality classes for charge averaged pions, kaons, and 



protons, as well for tt" [sJ] and (j) 52 1. We use previ- 
ously published PHENIX data on identified hadrons in 
p+p collisions [13] to evaluate the i?AA- The Raa data 
are limited in pT reach by the p+p data. As with the 
Rcp, the pions and kaons exhibit a suppression pattern 
in the -Raa- Additionally, a significant and monotonic 
centrality dependence is observed, with the suppression 
decreasing as the collisions become more peripheral. This 
is consistent with what is seen for neutral pions [13, H^j ■ 
The proton i?AA shows no suppression in the interme- 
diate Pt region and in fact reaches a maximum value 
above unity between 2-3 GeV/c. For p-r > 3 GeV/c, 
the proton i?AA values decrease and a suppression pat- 
tern emerges. The proton Rqp decreases more slowly 
than the central proton -Raa, which is simply because 
of the still considerable modification in the peripheral 
bins. The trend appears to be that the proton Rcp and 
-Raa decrease steadily while the pion Rc p and -Raa hold 
steady, suggesting that these values for pions and protons 
may eventually merge. The proton -Raa for the 0-10% 
centrality bin shown here exhibits reasonable qualitative 
agreement with the K + p Raa for 0-12% centrality re- 
ported by STAR [ig. 

While the centrality dependence of the i?AA for the 
pions and kaons is strong, it is quite weak for the pro- 
tons and the different centralities are consistent within 
the systematic uncertainties. This is consistent with the 
strong centrality dependence in the p/n ratios discussed 
above. The (p meson -Raa values are close to the val- 
ues for kaons, and significantly lower than the values for 
the protons, even though the 4> is much heavier than the 
kaon and it has roughly the same mass as the proton. 
This strongly suggests a baryon vs. meson dynamic, as 
opposed to a simple mass dependence, as would be the 
case for radial flow developed during the hadronic phase. 

The nuclear modification factor for d+Au collisions, 
RdA, is defined in a similar way as -Raa by 



RdA = 



(dN/dpT) 



d+Au+Au 



N^+^'^idN/dpTY+P' 



(6) 



Figure [12] shows RdA as a function of pT in different cen- 
trality classes for charged averaged particles. As with the 
-Raa, we use previously published PHENIX data on iden- 
tified hadrons in p+p collisions [i^l to evaluate the RdA- 
The p+p data limit the pt reach of the RdA- The charged 
pion exhibits a small modification above pt of 1.0 GeV/c 
and is consistent with nonmodification within the sys- 
tematic uncertainties. This is consistent with previous 
measurements of neutral pions [55, 58\. The charged 
kaon agrees with the charged pion within the system- 
atic uncertainties. The (f> meson exhibits no apparent 
modification. 

On the other hand, the protons show a very large 
and strongly centrality dependent Cronin enhancement, 
reaching a factor of 2 in the most central collisions at 
intermediate pt- Even in the 40-60% centrality class the 
enhancement is a factor of 1.5. For the most periph- 
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eral bin the enhancement is much smaller, at a factor of of the proton R^a is in fact very similar to the significant 
about 1.1-1.2, and is close to unmodified, similar to the centrality dependence of the p/n ratio, and these two ob- 
other particle species. This strong centrality dependence servables are likely driven by the same mechanism. Also 
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apparent in the RdA is that the enhancement for protons 
begins to fall off at 3.0 GeV/c and steadily drops with in- 
creasing pt, appearing nearly unmodified at the highest 
Pt points. 

The RdA of TT, K , (f), and p show significant dependence 
on the number of valence quarks, and no dependence 
on particle mass. That the baryon RdA is quite differ- 
ent from that of the mesons suggests that recombination 
plays a role in particle production in d-l-Au collisions as 
well as Au-|-Au. The kaon RdA is consistent with the pion 
RdA ) in contrast to Raa where the kaons are consistently 
above the pions. This is consistent with the K/tt ratio 
discussed above and indicates that there is no discernible 
strangeness enhancement within uncertainties in d-l-Au 
collisions. 



D. Comparison of Peripheral Au+Au to Central 
d+Au 

Motivated by the remarkable similarities between pe- 
ripheral Au-|-Au and central d-l-Au collisions, we now 
compare the two directly. Figure [T51 shows the K/tt ratio 
and Fig. ll4l shows the p/ir ratio in peripheral Au-|-Au and 
central d-l-Au collisions plotted together. In both cases 
the ratios are completely consistent with each other be- 
tween the different collision species, suggesting that the 
particle production mechanisms in peripheral Au4-Au 
and central d+Au collisions are quite similar. 
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FIG. 15; (Color online) Ratio of invariant yield of particles in 
peripheral Au-|-Au (60-92%) to central d+Au (0-20%) colli- 
sions as a function of pt- 

Figure [15] shows the ratio of the spectra in peripheral 
Au-|-Au to central d-fAu collisions for pions, kaons, and 
protons. Also plotted is the ratio for neutral pions, as 
determined from the data in Refs. [53,[5^. The ratio for 
TT^ shows excellent agreement with the tt^ ratios. 

We note that A^coii values are within 2% of each other 
and that the Apart values are within 4% of each other. 
Both of these differences are well within the associated 
systematic uncertainties of the Glauber MC values. No 
scaling is applied, but a scaling by the respective A'coii or 



A^part values would change the data imperceptibly. The 
ratios tend to the same value of roughly 0.65 for each 
particle species at and above 2.5-3 GeV/c. This univer- 
sal scaling is strongly suggestive of a common particle 
production mechanism between peripheral Au+Au and 
central d+Au collisions. It is also interesting to observe 
that at the lower px , where this ratio rises strongly, there 
is minimal mass or particle type dependence. 

Given the fact that both A"coii and A'part are essentially 
the same in these two systems, any quantity or physical 
effect that scales with either of these should be the same 
in each system, and thus should cancel almost exactly 
in this ratio. If we make the simple assumption that 
most or all of the cold nuclear matter effects scale with 
Acoii or A'part, then those effects are completely canceled 
in this ratio, leaving only the hot nuclear matter effects 
present in the peripheral Au+Au collisions. This could 
mean that this ratio being less than unity is attributable 
to the parton energy loss in peripheral Au+Au. This pic- 
ture is consistent with the findings in this paper and else- 
where that the Raa of mesons indicates parton energy 
loss in the medium even in peripheral Au+Au collisions. 
It is striking, then, that this ratio is independent of par- 
ticle species, which is suggestive of similar energy loss 
effects even for protons. This indicates that the baryon 
enhancement mechanism is the same in both systems. 

However, although A'coU and A'part are consistent for 
the two systems, there is an inherent participant asym- 
metry that needs to be taken into account. In the case of 
peripheral Au+Au collisions, one has a scenario in which 
7 or 8 nucleons on the edge of one Au nucleus collide 
against 7 or 8 nucleons on the edge of the other Au nu- 
cleus. On the other hand, in the case of central d+Au 
collisions, one has a scenario in which the 2 nucleons of 
the deuteron collide against 13 nucleons in the center of 
the Au nucleus. This introduces several additional fac- 
tors that need to be considered. For example, the partic- 
ipant asymmetry produces a rapidity shift in the particle 
production 59] . This may explain a deficit of soft parti- 
cles at low Pt in d+Au collisions at midrapidity, which 
in turn would explain why the ratio trends up at low px ■ 
We also note this low px region where the ratio rises is 
where hydrodynamics effects are known to be important 
in Au+Au collisions. It is possible that there are collec- 
tive fiow effects in d+Au collisions as well, as suggested 
by the recent results reported in Refs. f26l - l29| . A full 
viscous hydrodynamics model comparison is warranted. 

Another issue to consider is the modification of par- 
ton distribution functions (PDFs) in nuclei. These nu- 
clear PDFs (nPDFs) are known to be modified from the 
PDFs of single nucleons 0, HOj ■ The experimentally mea- 
sured nPDFs are averaged over the entire nucleus and 
are typically compared to the deuteron PDF to deter- 
mine the nuclear modification. The binary collisions in 
peripheral Au+Au involve two nucleons which have mod- 
ified nPDFs. On the other hand, the binary collisions in 
central d+Au involve an approximately unmodified nu- 
cleon from the deuteron and a modified nucleon from the 
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Au nucleus. Physical obscrvables sensitive to the nPDFs 
would then be expected to be different for the two sys- 
tems. However, it is possible that the nucleons in the 
more diffuse outer region of the nucleus have a differ- 
ent modification from those in the denser center. There- 
fore, it is not possible to make any model-independent 
quantitative statements about the differences between 
the nPDFs in these two systems. 

V. SUMMARY 

In summary we present a highly detailed and system- 
atic study of identified charged hadron spectra and ratios 
as a function of pT and centrality for Au+Au and rf+Au 
collisions at y^s^ 200 GeV. As has been reported pre- 
viously, we find a baryon enhancement present in both 
systems. In d-|-Au collisions, the Cronin enhancement 
has long been known to be stronger for baryons than 
for mesons. However, for the first time a study with 
enough statistical and systematic precision presents clear 
evidence for a strong centrality dependence of this ef- 
fect. In Au+Au collisions the baryon enhancement has 
been attributed to parton recombination as the mode of 
hadronization. A version of the recombination model has 
been applied to ri+Au collisions as well, which repro- 
duces the baryon vs. meson differences. The present data 
strongly suggest that further theoretical investigation is 
warranted. Given the excellent statistical precision of the 
present data set, a direct comparison between the two is 
made for the first time. Specifically, a ratio of the spectra 
in the most peripheral Au-|-Au and most central d-l-Au 
collisions is measured. These two systems have nearly 
identical values of both A^coii and iVpart- Therefore, a 
direct comparison between the two cancels out a large 
number of physical effects. We conclude that the baryon 
enhancement present in both systems is likely driven by 
a common hadronization mechanism. 
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